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ABSTRACT We study a physical model for the formation of bud- 
like invaginations on fluid membranes under tension, and apply 
this model to caveolae formation. We demonstrate that budding 
can be driven by membrane-bound inclusions (proteins) provided 
that they exert asymmetric forces on the membrane that give rise 
to bending moments. In pctrticulcir, Caveolae formation may not 
necessarily require forces to be applied by the cytoskeleton. Our 
theoretical model is able to explain several features observed ex- 
perimentally in caveolae, where proteins in the caveolin family are 
known to play a crucial role in the formation of caveolae buds. 
These include (i) the formation of caveolae buds with sizes in the 
lOOnm range (ii) that a fairly large variation of bud shape is ex- 
pected (iii) that certain N and C termini deletion mutants result 
in vesicles that are an order of magnitude larger. Finally, we dis- 
cuss the possible origin of the morphological striations that eire 
observed on the surfaces of the caveolae. 



1 Introduction 

It has long been understood that invaginations form spontaneously on cell 
membranes [1]. These invaginations, which eventually separate from the 
membrane as mature, membrane bound vesicles, play an essential role in 
cellular trafficking and signalling[2, 3]. The mechanism by which such in- 
vagination is controlled is still far from fully understood, although it is now 
widely accepted that certain membrane-bound proteins, including clathrin 
and caveolin, play an important role. The formation of clathrin coated pits is 
thought to be driven by the controlled geometric aggregation of clathrin into 
rigid scaffolding which force the membrane to curve [4, 5]. The mechanism 
for formation of the second most common class of membrane invaginations, 
known as caveolae, is less well understood. Caveolae, which are less morpho- 
logically distinct than clathrin coated pits, resemble Q-shaped invaginations 
with a typical size of the order of 100nm[6, 7, 8] . They are present at high con- 
centrations on primary adipocytes, fibroblasts, muscle cells, pulmonary type 
1 cells as well as endothelial cells, and perform a variety of function ranging 
from signal transduction to intracellular transport[9, 10]. A "striated coat" 
can be seen on the cytoplasmic side of the caveolae membrane. It is believed 
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to reflect the organization of a recently discovered class of membrane-bound 
proteins, called caveolins, which are crucial to the formation of cavcolac[3]. 

The protein caveolin has an hairpin structure, with a short membrane 
spanning sequence, flanked by two hydrophilic termini, both found on the 
cytoplasmic side of the cell membrane: a 101 amino-acid polypeptide N- 
terminus tail (polymer), and a shorter (44 a-a) C-terminal, which is strongly 
attached to the membrane[10]. These caveolin molecules are typically found 
in small aggregates of 15-17 moleculcs[7, 11], the aggregation being driven 
by residues of the N-terminal located close to the membrane. Furthermore it 
is believed[10] that there exist some specific C-terminal - C-terminal attrac- 
tions, which are responsible for the organization of the protein aggregates 
at the surface of the caveolae membrane. Mutational analysis of caveolin- 
induced vesicle formation have been recently performed[12] and are discussed 
in relation with our theory in section 7 below. 

Caveolae are now thought to influence cell physiology in many ways, 
including growth and cell division, adhesion, and hormonal response [13]. 
These invaginations have been associated with the formation of lipid rafts[14], 
glycosphingolipid- and cholesterol-enriched microdomains within the plasma 
membrane of eucaryotic cells. Their ability to perform many different tasks 
might be achieved by their involvement in reporting change in membrane 
composition by signal transduction to the nucleus. It may also be connected 
to their regulation of signal traffic in response to extracellular stimuli, in- 
cluding mechanical stress [15]. 

From a physical point of view, spontaneous vesicle formation has been ob- 
served in vitro by adding amphiphihc polymers to various lipid systems[16]. 
It can be viewed as an example of the so-called curvature instability of fluid 
membranes containing inclusions, predicted to occur for inclusions that lo- 
cally influence the membrane curvature [17]. There have been physical studies 
of the inclusion-induced budding of vesicles[18, 19, 20] and works on the effect 
of single[21] and distributed [22] polymers grafted on membranes. 

Our aim is to study the effect of small inclusions, such as proteins, that 
affect the shape of the cell membrane. We assume that this "foreign" ob- 
ject exerts a force on the membrane, which may be due either to entropic 
effects, similar in origin to the pressure exerted by a gas onto the walls of 
its container, or to specific mechano-chemical forces. Throughout we will 
attempt to compare our rather general theory with the specific phenomenon 
of caveolin-mediated formation of caveolae. The fact that membrane-bound 
objects exert a force on the membrane arises naturally from theories that de- 
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scribe polymers grafted to surfaces. These have been extensively developed 
over the last decade or so, based on early ideas due to de Gennes[23] and 
others. The forces exerted by membrane-bound inclusions, as well as there 
interactions, have been calculated in certain ideal situations (e.g. idealized 
polymers on a tensionless membrane) [24, 25]. In what follows, we will ana- 
lyzed arbitrary force distributions, which allow for the description of specific 
inclusions, such as the caveolin aggregates. We also focus on tension bear- 
ing membranes, a situation that we believe more closely approximates the 
plasma membrane of the cell. 

The deformation of a membrane subjected to an arbitrary force distribu- 
tion, and the subsequent membrane-mediated interaction between two such 
distributions, are computed in sec. 2. The description of possible physical 
models for the force distribution follows in sec. 3. A physical theory for the 
several levels of protein self-organization at the cell membrane is presented 
in section 4, followed by our results for different force distributions (sec. 5). 
We then briefly comment on possible physical mechanisms for the peculiar 
phase behavior (stripe formation) at the bud surface (sec. 6), and conclude 
in section 7. 

2 Membrane response to an arbitrary force 
distribution 

The deformation energy of a membrane involves its surface tension 7 and 
bending rigidity k. Cells commonly adjust their surface tension to a set 
value via a mechanism known as Surf ace- Area- Regulation [26]. Hence mem- 
brane phenomena over sufficiently long timescales effectively occur at con- 
stant surface tension. It is also known that the composition of biological 
membranes exhibits spatial variations. Caveolar membranes for instance, 
show a high cholesterol content [13], the precise biochemical role of which is 
not yet entirely clear. Prom a physical point of view, it is known that choles- 
terol increases the local rigidity of the membrane[27, 28]. Local variations of 
membrane rigidity are not included in the following model, but some (lim- 
ited) information on the impact of cholesterol on caveolae at the physical 
level can be obtained by examining the effect on uniform changes in k across 
the whole membrane. 

Initially, we restrict our analysis to a membrane that is weakly deformed 
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by the presence of the inclusions. We proceed by writing down the free energy 
of an infinite planar fluid membrane as a standard expansion in powers and 
gradients of the membrane displacement u{r) from its flat (unperturbed) 
position[29]. This is modified to include the leading order term arising from 
the coupling to u of the applied pressure distribution /(r), which will arise 
from the action of inclusion (s) [30]: 



(1) 



Typically values for phospholipid bilayers are[27] k ~ 20A;bT*. The sur- 
face tension is reported to be in the rangc[31] 7 k, 10~^ — 10~^ pN/nm. 
The interplay of surface tension and bending rigidity defines a characteristic 
lengthscale k"^ = \/ nj^ ~ 30 — 90nm. 

Minimization of this energy results in the equilibrium membrane dis- 
placement M(r), and is reported in more detail elsewhere [32]. We find ■u(r) = 
J G{x — r')/(r')(i^r' with the Green's function (the response to a point force) 
given by 

i^o(fc|r — r'l) + log/c|r — r 



G(r - r') = 



(2) 



2777 

where Kq is a modified Bessel function [33], which decreases exponentially 
over a size . The membrane displacement is discussed further in section 
3. 

Overlap of displacements due to neighboring inclusions lead to membrane- 
mediated interactions between them. The interaction potential $(r) between 
two similar inclusions separated by a vector r is obtained by inserting the 
total force distribution /(r') + fir' + r) into Eq.(l) and identifying the r- 
dependence of the resulting energy (see [32] for the general theory). The 
interaction energy reduces to 

$(r) = - y" d\' J dh"f{r')f{r")G{r - r' + r") (3) 

where G{r) is the real space Green's function given by Eq.(2). 

If the inclusions have a circular symmetry (/(r') = f{r')) and do not 
overlap (r > 26 where b is the spatial extent of the force) , wc have been able 

*It is usually helpful to compare energies to the energy available from thermal fluctu- 
ations ksT, where ks is the Boltzmann constant and T the temperature 
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to determine the interaction potential exactly in an analytic form 



*(r) = ^CMkr) (4) 

where C = 2Tir'dr'ip{r-')lQ{kr') characterizes the strength of the interac- 
tion (Jo is another modified Bessel function [3 3]). The interaction is every- 
where repulsive in the regime of validity r > 26. 



3 Models for membrane-bound proteins 

Up to this point, we have been able to avoid making any but a few rather 
general assumptions about the form of the force exerted by the membrane 
inclusions. We will now proceed to consider some specific models for the 
force distribution. We do this both to make possible the later quantitative 
comparison with experiments and to demonstrate how such forces are ex- 
pected to arise on general physical grounds. The force distributions and the 
subsequent membrane deformations (from Eq.(2)) are shown in Fig.l. 

3.1 Flexible polymers 

In this section we treat caveolin proteins as flexible, linear polymer chains, 
anchored to the membrane. The idealized picture enable us to extract an 
analytic estimate of the force distribution. In the language of polymer 
physics [2 3], the caveolin homo-oligomer can be viewed as a brush of Q ~ 16 
polymer chains, grafted by one end to a small patch of membrane of radius 
a. The flexible chains on average arrange themselves radially to form a hemi- 
sphere of radius h (Fig. la). Thus for radial distances a < r < b one finds a 
corona of randomly coiled polymer chains with a chain density that is larger 
near the core and smallest on the outskirts of the distribution. 

A central concept in the theory of polymer physics is the existence of 
a correlation length or 'blob size' ^(r)[23], which is roughly the distance 
between interchain contacts in the corona of the caveolin brush (viewed as 
a semi-dilute polymer solution). Each chain can then be thought of as a 
string of correlation blobs extending radially outwards, with small values of 
^ corresponding to large densities of monomers. The classical Daoud-Cotton 
model [34] takes advantage of the fact that the surface area of a hemisphere 
of radius r is approximately filled by Q close-packed blobs, to deduce the 
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scaling of the correlation length: ,^(r) = r^j2T{/Q. It can be shown [23] that 
the work done in generating each blob is ksT ^ independent of the blob size. 
Thus we may write the pressure in this region as the energy per blob divided 
by the volume of a blob 



A result that is valid for a < r < 6, and that is consistent with more detailed 
calculations [24, 25]. The physical origin of this pressure can be understood 
as simply due to the exchange of momentum due to collisions of the polymer 
chains with the membrane. It is therefore purely entropic in origin, as evi- 
denced by the overall ksT energy scale. For r > h the pressure is zero since 
the chains have finite length (in fact, it is exponentially small[24, 25]). The 
pressure in the core binding region is assumed constant and must involve a 
total force equal and opposite to that applied by the corona[30] 



.1 ifO<r<a 

/(^) = i if -<r<h (6) 

' if r>6 



Note that the strength of the force is characterize by the energy scale Eq 
f[r = 0)a^. For the Q-chain oligomer, it is Eq = 2kBT{Q/2nf/^{l - a/b) 
AkeT. 



3.2 Block Distribution 

We believe that the polymer brush model captures some of the fundamental 
properties of a collection of large hydrophilic proteins anchored to a biomem- 
brane, namely i) a. downward pressure exerted by the cytosolic portion of the 
proteins, combined with ii) an upward pull from the anchors (the hydropho- 
bic region of the proteins). It however employs rather strong assumptions 
(random coil configuration, absence of internal structure, and large size of 
the polymer chains) which are certainly not satisfied for the protein caveolin. 
The simplest example of a general force distribution that satisfies the criteria 
above is a "block distribution" , for which the force exerted by the hydropho- 
bic anchors (between < r < a) and the hydrophilic sections {a < r < b) 
are both constant: E^/a? = fr^a = {b'^/a^ — ^)fa<r<b- The membrane de- 
formation for such distribution is larger than for the brush distribution for 
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the same strength, as characterized by the energy Eq. since the force is not 
concentrated near the center of the distribution (Fig.lc). 



4 bud structure and morphology 



Caveolae formation involves a hierarchy of self organization, ranging from the 
nanometric scale (oligomers of ~ 15 particles and size ~ 5 nm) up to buds of 
radius ~ 100 nm. We review briefly the theoretical framework of thermal self- 
organization [2 9], and give insights on the caveolin homo-oligomer formation, 
which we view essentially as a micellization in two dimensions. We then 
describe in some details the formation of caveolae buds. 

Consider a solution of particles of average surface fraction 0, that can 
exist either as isolated entities or in larger aggregates (homo-oligomers) of p 
particles and of energy Fp = pfp. The concentration Cp of p-sized aggregates 
follows a Boltzmann law[29]: Cp ~ q-{Pp-w)/i^bT ^ where /i is the chemical po- 
tential of the particles, usually fixed by the average concentration 0. There 
is an energetic tendency to form aggregates if the energy per particle fp is 
(at least in some regime) a decreasing function of the aggregation number p. 
It overcomes the entropic dispersive effect beyond a critical value of (the 
critical aggregation concentration, or c-a-c), usually defined as the concen- 
tration at which the density of aggregates is equal to the density of isolated 
particles. At the c-a-c and above the average size p* of the aggregates is the 
one that minimizes the energy per protein in the aggregate. The root mean 
squared deviation Ap from the average depends upon the steepness of the 
energy variation around that minimum. Expressed in the form of equation, 
these conditions yield (for p* ^ 1): 



The driving force for homo-oligomcrization is an attractive interaction 
between specific motives on the N-terminal of the protein[ll]. We proceed 
by assuming that all proteins in the interior of the oligomer experience a mu- 
tual attraction, and contribute to the oligomer energy Fp via a negative linear 
term —fx'p. Proteins in the outskirts of the oligomer on average experience 
less attraction, as they have less neighbors. They increase the energy Fp by a 
factor +7r/3y^, where /? is the energy loss {per protein) for being on the out- 




dp 
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skirts^ Bringing proteins together also leads to steric and/or entropic repul- 
sion (crowding). For a polymer brush for instance, the latter contribution to 
Fp is +ap^^'^ (see Eq.(5)). Solving Eq.(7) with Fp/p = fp = a^+vr/?/^— /x', 
the experimental observation (oligomers containing 14— 16 proteins[ll]) are 
consistent with a — 2kBT and /3 = lO/csT. Both numbers are physically 
reasonable: the frustration a which results in bringing many proteins close 
to one another might be expected to be of order ksT per protein, and the cal- 
culated energy loss f3 for proteins on the outskirts of the oligomers is of order 
an hydrogen bond energy {10k bT) per proteins. Although quite crude, this 
model provides a thermodynamic description for the first level of self organi- 
zation in caveolar membranes, which is able to reproduce the experimental 
observation on cavcolin oligomerization, namely p* ~ 15 and Ap ~ 2. 

The formation of the caveolae themselves can be described more accu- 
rately. As it involves a large number of oligomers, a precise description at 
the molecular level is less crucial. We model the Q-shaped invagination by 
a closed sphere of radius R. The bending moments exerted by the protein 
oligomers, which are all on the cytoplasmic side of the membrane, drives 
the bud formation, expected to occur above a critical budding concentration 
(c-b-c) of oligomers (see Fig. 2). For small concentration < (pcbc (Fig-2a), 
the membrane is uniformly covered by oligomers and remains almost fiat. 
Buds start forming as the concentration increases, and outnumber isolated 
oligomers at the c-b-c (Fig. 2b). If the concentration is increased further 
(Fig. 2c), then the concentration of isolated oligomers, and the bud size, re- 
main almost constant (0i = (pcbc), while the number of buds increases. 

The free energy per membrane inclusion in the bud ft, (Eq.(8) below) 
contains several contributions. Energy is gained if the membrane curves 
to accommodate the deformation imprinted by the cavcolin oligomer. A 
membrane curving away from the caveolin aggregate is favored (first term, 
RHS of Eq.(8)). In the limit of small curvature, the energy reduction per 
oligomer is of order —D^r/R (Eq.(9)). On the other hand, bud formation 
costs an energy that depends on the bending rigidity and surface tension 
(second and third terms, RHS of Eq.(8)), and leads to higher local oligomer 
concentration, modifying the pair interaction energy (fourth term, RHS of 
Eq.(8)). This interaction is characterized by the second virial coefficient B2 
(Eq.(9)), and involves the full interaction potential V{r) (see sec. 5). The 
last term in the RHS of Eq.(8) is the mixing entropy of a gas of membrane 

Hhe equivalent of the surface tension of a liquid 



8 



inclusions on a lattice. 



Drr , ^ KSi 7Si ,kBTB2 

+kBT (^log + - 1) log(l - (8) 

where si = 7r6^ is the oligomer area, and with 

Drr = J ^rV(r) B2 = j — [l-e^BT ^ (9) 

Minimizing this energy with respect to R gives the optimal bud radius 
R* = Rmin/4>, where Rmin = 4:KSi/Drr corrcspouds to the minimal radius of 
a bud densely packed with caveolin. Further energy minimization leads to 
the optimal amount of protein 0* recruited in the bud, defined by: 

log (1 - (j)*) 7Si 2/s:gi 



mvn 



This equation has a clear physical meaning. The protein coupling to the 
membrane curvature effectively reduces the second virial coefficient by an 
amount B^^^ = 2ks\/ {ksTR^^^, which indicates an attraction between 
oligomers [17]. 

The optimal concentration of Eq.(lO) corresponds to an energy minimum 
/j*. At the c-b-c (Eq.(7)), it is equal to the oligomer chemical potential: 
fb — l^cbc- The latter can be related to the concentration 0i of isolated 
oligomer via ii = log(0i/(l — 0i)) + 25201, which is the chemical potential 
of a gas on a lattice with pair interaction. The equation defining the critical 
budding concentration is: 



log 7 h 2 (f)cbc = 

1 - (Pcbc Si 

log- -- + 2 (11) 

1 — 0* Si 

The mean variation of the radius AR^ = {R^) — R*'^ can be approximately 
calculated by using a steepest descent method to calculate moments of the 
bud size distribution (Eq.(7)). We find {AR/R*f = ksT / {IQt:k)B'J {B'^ - 
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B2-'^), where 3'^ = B2 + Si/(20*(1 - 0*)) > includes both the interaction 
between brushlets and the entropic contribution-'-. The mean radius varia- 
tion shows the signature of the membrane curvature instabihty mentioned 
earher[17]. If the couphng between membrane and inclusion is sufficiently 
strong: B^^^ > B'2, or D"^^ > 8kBTKB'2, then small fluctuations of any 
lengthscale are unstable and the mean variation of radius becomes large. 
The actual dispersion in bud size depends on how close we are to the insta- 
bility. It is of order 6% for the parameters used below. Note that variations 
in shape that conserve the mean curvature of the membrane should be larger, 
as they only cost a fraction of the energy penalty corresponding to variation 
of the bud global size. Prom electron micrographs of caveolar membranes, 
the projected radius variation is of order 20%[6] (see also Fig. 13.48 in [1]). 

Numerical calculation of the bud radius and protein concentration in the 
membrane is shown in Fig. 3 upon variation of the surface tension for different 
values of the coupling strength (for an attractive energy Eatt — 0.5 /cbT", see 
sec.5). A strong variation in bud size is observed for small surface tension. 
At larger tension, the radius is almost insensitive to 7. Bud formation is 
however less favorable, as can be seen from the increase of the c-b-c. Our 
model also predicts the existence of a critical point [35], hence a possible 
coexistence of buds of different radius, connected to the curvature instability 
studied by Leibler[17]. We will not discuss this further here, as it is probably 
not relevant to the problem of caveolae formation. 

We have derived the bud morphology as a function of two variables which 
depend on the actual shape of the force distribution: B2 and Drr- In order 
to make quantitative prediction, we study below two "extreme" force distri- 
butions. 

5 Results for various force distributions 

In this section we discuss the results above for the polymer brush model 
and the block model. The force distributions involve 3 parameters. The 
lengthscales a and b can be measured experimentally: a ~ 2nm and b ~ 
5nm. The energy scale of the force Eq can be calculated for the brush 
model, and will be estimated for the block model(Eq.(6)). We present results 
for the minimal radius Rmin (connected to the force moment Drr) and the 

■tif B2 < 0, the inclusions spontaneously demix on the flat membrane 
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membrane mediated interaction $(r). The excluded volume B2 involves the 
full interaction potential and is discussed at the end of this section. 

The calculation results in a unified description of the force distribution 
via the energy scale Eq and the ratio of size of the protein aggregate over 
deformation range ka (for ka <^ 1). The bud radius is Rmin/b — ckk/Eq with 

16a 16 / a\ , . 

Oiblock = ~ 6 abrush ^ "3" + j ^ '^^^f 

The oligomer-oligomer interaction is given by Eq.(4), and involves the force 
moment ( = —(3{kayEQ/a with 

= ^ = 2.5 PbrusH = ^^^^l = ^^^^ 

Brush distribution 

For ideal gaussian chains, the energy scale given by Eq.(6) is Eq ~ 
((5/27r)^/^ ~ AksT. Most of the force is concentrated near the center of the 
distribution, and has a small effect on the membrane. Fig.l shows that the 
deformation is quite small (of order OAnm). However, collective effects lead 
to the formation of fairly small buds of minimum bud radius Rmin — 300nm 
(much larger, however, than the caveolae). The membrane mediated inter- 
action between protein aggregates is very small, of order 10~^A;bT. 

Block distribution 

The block distribution is probably more relevant to the case of stiff, short 
proteins such as cavcolin. We choose the strength of the force so that each 
protein contribute to of order ksT of energy: Eq ~ 10 — 15k bT, which im- 
poses a displacement u{r = 0) = 2nm (Fig.l). The corresponding minimum 
bud radius Rmin — 60nm, is comparable to the radius of caveolae. Radius 
variation with surface tension is shown in Fig.3. The interaction potential 
is $(r) = 0.02KQ{kr) ksT. We believe that, although small, this repulsive 
interaction might be responsible for the remarkable phase behavior of the 
proteins at the surface of the buds (sec. 6). 

Second virial coefficient and critical budding concentration 
The virial coefficient B2 defined by Eq.(9) involves the full interaction 
potential V{r), including hard core repulsion and the membrane mediated 
physical interaction Eq.(4). Moreover, we have experimental evidence[10, 36] 
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that there exists short range specific attractions between the protein side 
chains (C-termini). To describe this attraction, we adopt the exponentially 
short range form Vatt = —Eatt&~^^''^^^^i where Eatt ~ ksT is the strength of 
the attraction, of range the size of the oligomer h. The short range attraction 
acts to increase the oligomer density inside the invagination. The resulting 
buds are crowded with proteins ~ 0.8, and are quite small, with a radius 
of order R ~ 70nm - see Fig. 3. 

6 Microphase separation at the bud surface 

One peculiar feature of the caveolae is their striated texture, believed to cor- 
responds to alignment of protein ohgomers at the surface of these buds [6] (see 
also Fig. 13.48 in [1]). This finding is particularly striking as it is not trivial 
to understand how radially symmetrical oligomers may organize themselves 
into non-symmetrical phases. We argue that the stripe phase might be a 
signature of the membrane mediated repulsion between protein aggregates. 
Molecular dissection of the caveolin protein has shown that the oligomers 
interact attractively via the third distal region of their C-termini [36]. This 
attraction may lead to gas/liquid phase separation of the caveolin oligomers, 
which results in dense membrane patches (the liquid) coexisting with less 
dense regions (the gas). Our situation is more complex, as we have shown 
the existence of an additional, membrane-mediated, longer range repulsion 
between oligomers. It has recently been argued at the light of computer 
simulation [37] that under this long range repulsion, the gas and liquid phases 
are broken into microdomains (circles at small concentration, and stripes 
for higher concentration). This is because large aggregates are costly due 
of the long range repulsion, while small aggregates (circles at low density, 
and stripes for larger density) are favored by the short range attraction. A 
linear stability analysis for a solution of particle interacting with an expo- 
nentially short range attraction (of strength Ea and range the particle radius 
6), added to the membrane mediated repulsion of Eq.(4) (of strength E'^), 
shows[32] that periodic arrays of dense and dilute regions are expected for 
strong enough repulsion E^ > ^Ea{khY. Because of the long range of the 
repulsion (for oligomers, h ~ 5nm, ~ 50nm), a repulsive interaction as 
low as 10~^A;bT between protein oligomers can indeed produce a well or- 
dered phase for an attractive interaction of order ksT. The structures have 
a typical size 27r/g* which is, to lowest order, independent of the range of 
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the repulsion: q*b = {2Er/2>EaY/^ — \- defines dense regions of 

size of order a few oligomer diameters, which compares quite well with the 
experimental observations. 

7 Conclusions 

We have constructed a theory for the formation of bud-hke invaginations 
on fluid membranes, which we have compared with experimental data for 
caveolae found on the surfaces of cells. Our description incorporates the 
effect of both membrane rigidity and tension. 

Proteins in the caveolin family are known to play a crucial role in the 
formation of caveola, by forming homo-oligomers which concentrates in the 
buds. We argue that asymmetrically anchored membrane proteins (or protein 
oligomers for caveolin) can apply forces to the membrane. We examined 
several models for the origin and magnitude of these forces, which may be 
purely entropic in origin or may result from stronger interactions. Such forces 
act to exert bending moments on the membrane and drive the formation of 
bud-like structures, for which we are able to make theoretical predictions. 
Our model correctly reproduces the size of the buds (lOOnm). and provides 
a physical explanation for the origin of the morphological striations observed 
on their surface. 

Our results also shed hght on several experimental observations concern- 
ing the function of caveolae and the result of caveolin mutation. It has 
recently been suggested that caveolae-like domains play a critical role in 
the mechanosensing and/or mechanosignal transduction of the extracellular 
signal-regulated kinase pathway[15]. We predict that while the membrane 
tension has little effect on the size of the buds (Caveolae indeed have a sim- 
ilar structure in different kind of cell, which possibly bear different tension), 
the amount of caveolin protein required to observe bud formation does in- 
crease strongly with surface tension. Increasing the tension (either via a 
shear stress or by direct cell manipulation) may result in the disappearance 
of the buds if the available protein amount is insufficient. This is a testable 
prediction. 

We have also described how the control of the morphology of caveolae 
can be achieved in a number of ways, including by the level of cholesterol 
in a membrane. High cholesterol content in membrane is known to result in 
higher membrane rigidity, driving an increase in the bud radius. 
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Finally, caveolae formation in mutant caveolin systems[12] also provides 
verification of our theory. Mutants which lack the self-attractive segment 
of the N-termini (responsible for the formation of the homo-oligomers) are 
still competent to drive vesicle formation, but result in much larger buds 
R ~ liim. This is consistent with the fact that the force exerted by isolated 
proteins should be about 10 times smaller than the force exerted by oligomers 
resulting in a ten-fold increase of the bud radius. Mutants which lack the 
mutually attractive C-terminus result in similarly larger buds. Within our 
theory, this mutation results in a weaker oligomer-oligomer attraction, hence 
in a lower density of caveolin in caveolae and therefore larger buds. However, 
our theory predicts that ohgomer attraction should not strongly influence 
the bud size. A natural conclusion would be that the C-termini contribute, 
either directly or indirectly, to the force exerted on the membrane, and that 
this force is reduced in mutants. 
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Figure Captions 



FIG. 1. (a) Sketch of the blob model for the anchored protein aggregate and 
(b) force distribution for the two models used in the paper: brush distribution 
(dashed, corresponding to (a)) and block distribution (sohd). The membrane 
is pushed down by the "corona" of grafted polymers out to r — b — 5nm and 
is pulled upwards by the anchored "core" inside r' = a = 2nm. (c) The 
corresponding membrane deformation u{r) in unit -^a for = 3Qnm. The 
brush distribution has a weaker effect on the membrane because the force is 
mostly concentrated near its center (r — 0). For aggregates residing on the 
cytoplasmic face of the membrane, including caveolin homo-oligomers the 
cell interior would be above the membrane. 

FIG. 2. Sketch of bud formation upon increase of oligomer concentration, 
(a) Below the critical budding concentration (c-b-c), the membrane is uni- 
formly covered by isolated oligomers, (b) At the c-b-c, buds have formed and 
outnumber isolated inclusions, (c) Above the c-b-c, the size and shape of the 
buds remains the same, and their number increases with the concentration. 

FIG. 3. Variation of the caveolae preferred radius R* (in nm) with the surface 
tension for 2 values of the coupling strength Eq for a short range attraction 
of Eatt = O.Sfc^T between brushlets (see sec. 5). Two inset curves show the 
variation of the bud composition 4>* (solid) and critical budding concentration 
(dashed) with for the same range of surface tension. The other parameters 
are a — 2nm, b — Bnm, k — 20/cbT 
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Figure 1 - Sens & Turner 
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Figure 2 - Sens & Turner 
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